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Abstract 

We extend the Standard Model by adding two gauge-singlet Z2-symmetric scalar fields that 
interact with visible matter only through the Higgs particle. One is a stable dark matter WIMP, 
and the other one undergoes a spontaneous breaking of the symmetry that opens new channels 
for the dark matter annihilation, hence lowering the mass of the WIMP. We study the effects 
of the observed dark matter relic abundance on the annihilation cross section and find that in 
most regions of the parameters space, light dark matter is viable. We also compare the elastic 
scattering cross-section of our dark matter candidate off nucleus with existing (CDMSII and 
XENONIOO) and projected (SuperCDMS and XENONIT) experimental exclusion bounds. We 
find that most of the allowed mass range for light dark matter will be probed by the projected 
sensitivity of XENONIT experiment. 
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I. INTRODUCTION 



Cosmology tells us that about 25% of the total mass density in the Universe is dark 
matter that cannot be accounted for by conventional baryons [1]. Alongside observation, 
intense theoretical efforts are made in order to elucidate the nature and properties of 
this unknown form of matter. In this context, electrically neutral and colorless weakly 
interacting massive particles (WIMPs) form an attractive scenario. Their broad properties 
are: masses in the range of one to a few hundred GeV, coupling constants in the milli-weak 
scale and lifetimes longer than the age of the Universe. 

Recent data from the direct-detection experiments DAMA/LIBRA ^ and CoGeNT 
3|, and the recent analysis of the data from the Fermi Gamma Ray Space Telescope 
J|, if interpreted as signal for dark matter, require light WIMPs in the range of 5 to 
10 GeV [sj. Also, galactic substructure requires still lighter dark matter masses 0, 0]- 
In this regard, it is useful to note in passing that the XENONIOO collaboration has 
provided serious constraints on the region of interest to DAMA/LIBRA and CoGeNT 
assuming a constant extrapolation of the liquid xenon scintillation response for nuclear 
recoils below 5 keV, a claim disputed in [9|. Most recently the CDMS collaboration has 
released the analysis of their low-energy threshold data [10| which seems to exclude the 
parameter space for dark matter interpretation of DAMA/LIBRA and CoGeNT results, 
assuming a standard halo dark matter model with an escape velocity v^sc = 544 km/s and 



neglecting the effect of ion channeling ll|]. However, with a highly anisotropic velocity 



distribution, it may be possible to reconcile the CoGeNT and DAMA/LIBRA results with 



the current exclusion limits from CDMS and XENON 12] (se e also comments on p. 6 in 



13| about the possibility of shifting the exclusion contour in [lO| above the CoGeNT signal 



region). In addition, CRESST, another direct detection experiment at Gran Sasso, which 
uses CaW04 as target material, reported in talks at the IDM 2010 and WONDER 2010 
workshops an excess of events in their oxygen band instead of tungsten band. If this signal 
is not due to neutron background a possible interpretation could be the elastic scattering 
of a light WIMP depositing a detectable recoil energy on the the lightest nuclei (oxygen) 
in the detector [l4J. While this result has to await confirmation from the CRESST 
collaboration, it is clear that it is important as well as interesting to study light dark 
matter. 

The most popular candidate for dark matter is the neutralino, a neutral i?-odd super- 
symmetric particle. Indeed, they are only produced or destroyed in pairs, thus rendering 



the lightest SUSY particle (LSP) stable 15|. In the minimal version of the supersym 



metric extension of the Standard Model, the neutralino is a linear combination of the 
fermionic partners of the neutral electroweak gauge bosons (gauginos) and the neutral 
Higgs bosons (higgsinos). They can annihilate through a t-channel sfermion exchange 
into standard model fermions, or via a t-channel chargino-mediated process into W^W~, 
or through an s-channel pseudoscalar Higgs exchange into fermion pairs and they can 
undergo elastic scattering with nuclei mainly through scalar Higgs exchange 16||. If these 
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neutralinos were light, they would then be overproduced in the early universe and, in the 
minimal model, would not have an elastic-scattering cross section large enough to account 
for the DAMA/LIBRA and CoGeNT results due to constraints from other experiments 
such as the LEP, Tevatron, and rare decays 17H19 . 



Therefore, with no clear clue yet as to what the internal structure of these WIMPs is, 
if any, a pedestrian approach can be attractive. In this logic, the simplest of models is 
to extend the Standard Model by adding a real scalar field, the dark matter, a Standard- 
Model gauge singlet that interacts with visible particles via the Higgs field only. To ensure 
stability, it is endowed with a discrete Z12 symmetry that does not spontaneously break. 
Such a model can be seen as a low-energy remnant of some higher-energy physics waiting 
to be understood. In this cosmological setting, such an extension has first been proposed 



in 20| and further studied in 2l|] where the unbroken Z2 symmetry is extended to a 
global U(l) symmetry. A more extensive exploration of the model and its implications 
was done in |22|, specific implications on Higgs detection and LHC physics discussed in 
23 1 and one- loop vacuum stability looked into and perturbativity bounds obtained in 



24| . The work of 25[ considers also this minimal extension and uses constraints from 
the experiments XENONIO [26| and CDMSII 27 1 to exclude dark matter masses smaller 
than 50, 70 and 75GeV for Higgs masses equal to 120, 200 and 350 GeV respectively. 

In order to allow for light dark matter, it is therefore necessary to go beyond the 
minimal one-real-scalar extension of the Standard Model. The natural next step is to add 
another real scalar field, endowed with a Z2 symmetry too, but one which is spontaneously 
broken so that new channels for dark matter annihilation are opened, increasing this way 
the annihilation cross-section, hence allowing smaller masses. This auxiliary field must 
also be a Standard Model gauge singlet. 

After this brief introductory motivation, we present the model in the next section. We 
perform the spontaneous breaking of the electroweak and the additional Z2 symmetries 
in the usual way. We clarify the physical modes as well as the physical parameters. 
There is mixing between the physical new scalar field and the Higgs, and this is one 
of the quantities parametrizing the subsequent physics. In section three, we impose the 
constraint from the known dark matter relic density on the dark-matter annihilation cross 
section and study its effects. Of course, as we will see, the parameter space is quite large, 
and so, it is not realistic to hope to cover all of it in one single work of acceptable size. 
Representative values have to be selected and the behavior of the model as well as its 
capabilities are described. Our main focus in this study is the mass range O.lGeV - 
lOOGeV and we find that the model is rich enough to bear dark matter in most of it, 
including the very light sector. In section four, we determine the total cross section cr^g^^ 
for non relativistic elastic scattering of dark matter off a nucleon target and compare it 
to the current direct-detection experimental bounds and projected sensitivity. For this, 
we choose the results of CDMSII and XENONIOO , and the projections of SuperCDMS 



28| and XENONIT 29|]. Here too we cannot cover all of the parameter space nor are we 



going to give a detailed account of the behavior of cxdet as a function of the dark matter 
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mass, but general patterns are mentioned. The last section is devoted to some concluding 
remarks. Note that as a rule, we have avoided in this first study narrowing the choice of 
parameters using particle phenomenology. Of course, such phenomenological constraints 
have to be addressed ultimately and this is left to a forthcoming investigation, contenting 
ourselves in the present work with a limited set of remarks mentioned in this last section. 
Finally, we have gathered in an appendix the partial results regarding the calculation of 
the dark matter annihilation cross section. 



II. A TWO-SINGLET MODEL FOR DARK MATTER 

We extend the Standard Model by adding two real, spinless and Z2-symmetric fields: 
the dark matter field 5*0 for which the Z2 symmetry is unbroken and an auxiliary field 
Xi for which it is spontaneously broken. Both fields are Standard Model gauge singlets 
and hence can interact with 'visible' particles only via the Higgs doublet H. This latter 
is taken in the unitary gauge such that if^ = l/-\/2(0 /i'), where h' is a real scalar. 
We assume all processes calculable in perturbation theory. The potential function that 
includes 5*0, h' and Xi writes as follows: 

u = ^sl - f^h- - fx? + |s5 + + |xt + ^slh- + 'Jfslxl + ^h-xl (2.1) 

where mg, yU^and ^\ and all the coupling constants are real positive numbers. In the 
Standard Model scenario, electroweak spontaneous symmetry breaking occurs for the 
Higgs field, which then oscillates around the vacuum expectation value v = 246GeV (3o| . 
The field xi will oscillate around the vacuum expectation value fi > 0. Both v and vi 
are related to the parameters of the theory by the two relations: 



^ ~ ^ Xvi - 36Af ' ~ ^ Xvi - 36A? ' ^^'^^ 



It is assumed that the self-coupling constants are sufficiently larger than the mutual ones. 

Writing h' = v + h and Xi = ^^i + >S'i, the potential function becomes, up to an irrelevant 
zero-field energy: 

U = f/quad + f^cub + f^quar, (2.3) 

where the mass-squared (quadratic) terms are gathered in t/quad? the cubic interactions in 
f/cub and the quartic ones in t/quar- The quadratic terms are given by: 

f/quad = Im'X + \Mlh^ + \mI~SI + MlJiS,, (2.4) 

where the mass-squared coefficients are related to the original parameters of the theory 
by the following relations: 

ml = ml + + fvl, Ml = -^^ + ^v^ + ^v^ 

Ml = + ^v^ + ^vl- Ml = X,vv,. (2.5) 



Replacing the vacuum expectation values v and Vi by their respective expressions fl2.2p 
will not add clarity. In this field basis, the mass-squared matrix is not diagonal: there is 
mixing between the fields h and Si. Denoting the physical mass-squared field eigenmodes 
by h and 5*1, we rewrite: 



TT ^ 2o2 , ^ 21,2,-'- 2o2 

where the physical fields are related to the mixed ones by a 2 x 2 rotation: 



(2.6) 



cos 9 sin 9 



h 



(2.7) 



sin 9 cos 9 J \Si^ 

Here 9 is the mixing angle, related to the original mass-squared parameters by the relation: 



tan 261 



Ml -Mr 

and the physical masses in (12. 6p by the two relations: 



(2.8) 



mi 



Ml + Ml + e {Ml - Ml) J {Ml -Mlf + 4Mf, 



Ml + Ml-e {Ml - Ml) J (Ml - Mlf + AMf, 



(2.9) 



where e is the sign function. 

Written now directly in terms of the physical fields, the cubic interaction terms are 
expressed as follows: 



A 



(3) 



^(3) 

Ucuh = —^S^h H — —SqSi 



A(3) 
6 



6^2 2 ^ 



(3) 



x(3) 

^1 U2i 



(3) 



(2.10) 



where the cubic physical coupling constants are related to the original parameters via the 
following relations: 



X(3) 
(3) 

Vol 



\qV cos 9 + riQiVi sin 9; 
rjoiVi cos 9 — XqV sin 9; 



A(^) = Af cos^6' + -Aisin26'(t;icos6' + t;sin6') +?7it;isin^6'; 



(3) 

Af) 
Af 



?7if 1 cos'^ 9 Ai sin 29 (v cos 9 — vi sin 9) — \v sin^ 9; 



2 

Ai^i cos^ 6^ + - sin 29 [(2Ai - A) cos 6^ - (2Ai - r/i) vi sin 9] - Xiv sin^ 9; 
Xiv cos^ 6^ - ^ sin 29 [(2Ai - r]i) vi cos 9 + (2Ai - \)v sin 9] + XiVi sin^ 9. 



(2.11) 
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Also, in terms of the physical fields, the quartic interactions are given by: 



n A(4) Ti(^) x(^) r,(^) a(^) 

'/O r-4 I ^ 7 4 , 71 c,4 , •^O q2i2 , '/Ol 02 02 , ^01 027 



u, 



quar 



X(4) x(4) XV 



(2.12) 



where the physical quartic coupling constants are written in terms of the original param- 
eters of the theory as follows: 



(4 

m 

(4 
^01 

(4 
01 



A 



(4; 



A 



(4; 



A 



(4 



A cos^ 9 + - Ai sin^ 29 + r]^ sin^ 9- 

r^i cos^ 9 + ^ Ai sin^ 2^ + A sin^ 9] 

Ao cos^ 9 + ?7oi sin^ 9; 
r]oi cos^ + Ao sin^ 9; 



iVoi - Ao) sin 29; 



- - [(3Ai - A) cos^ 9 - (3Ai - 771) sin^ 9] sin 29; 
= Ai cos^ 2^ - ^ (2Ai - 771 - A) sin^ 29; 

= o [(^1 ~ 3Ai) cos^ ^ - (A - 3Ai) sin^ 9] sin 29. 



(2.13) 



Finally, after spontaneous breaking of the electrowcak and Z2 symmetries, the part 
of the Standard Model lagrangian that is relevant to dark matter annihilation writes, in 
terms of the physical fields h and Si, as follows: 

UsM = Yl (^'^/W + A1/7S1/7) + x^lhw^w^^ + xf^SiW-w+^ 
f 

+aS/^ {z,f + xf^s, {z,f + xf^h-'w^w^^ + xflslw;:w^^ 

+XniwhSiW;W+^' + A^>2 ^z^f + xf}sl {Z^f + Xni^hSi {Z^f . (2.14) 

The quantities m/, rriyj and are the masses of the fermion /, the W and the Z gauge 
bosons respectively, and the above coupling constants are given by the following relations: 

A 
A 
A 

Alt) = ?4 cos^ 9; Ait^ = ?4 sin^ ^; Am. = sin 29. (2.15) 
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III. RELIC DENSITY, MUTUAL COUPLINGS AND PERTURBATIVITY 



The original theory fl2.ip has nine parameters: three mass parameters (mo, /i, /ii), 
three self-couphng constants (?7o,A,?7i) and three mutual coupling constants (Aq, //oi, ^i). 
Perturbativity is assumed, hence all these original coupling constants are small. The 
dark-matter self-coupling constant tjq does not enter in the calculations of the lowest- 



order processes of this work [31[, so effectively, we are left with eight parameters. The 
spontaneous breaking of the electroweak and Z2 symmetries for the Higgs and Xi fields 
respectively introduces the two vacuum expectation values v and vi given to lowest order 
in (12.21) . The value of v is fixed experimentally to be 246GeV and for the present work, 
we fix the value of Vi at the order of the electroweak scale, say lOOGeV. Hence we are left 
with six parameters. Four of these are chosen to be the three physical masses mo (dark 
matter), mi {Si field) and nih (Higgs), plus the mixing angle 9 between Si and h. We 
will fix the Higgs mass to rrih = 138GeV and give, in this section, the mixing angle 6 the 
two values 10° (small) and 40° (larger). The two last parameters we choose are the two 
physical mutual coupling constants Ag'^'' (dark matter - Higgs) and rjj^^ (dark matter - Si 
particle), see (12.121) . 

In the framework of the thermal dynamics of the Universe within the standard cos- 
mological model j32|, the WIMP relic density is related to its annihilation rate by the 
familiar relations: 



flnh 



2 _ 1.07 X lO^Xf 



g*mpi (fl2CXann) GcV ' 



, 0.038mpimo (t;i2crami) . x 

Xf ~ m — = . [6.1) 

The notation is as follows: the quantity h is the Hubble constant in units of 
100km/(sxMpc), m-pi = 1.22 x lO^^GeV the Planck mass, mo the dark matter mass, 
Xf = mo/Tf the ratio of the dark matter mass to the freeze-out temperature Tf and (7* 
the number of relativistic degrees of freedom with a mass less than Tf. The quantity 
(^i20"ann) IS the thermally averaged annihilation cross-section of a pair of two dark matter 
particles multiplied by their relative speed in the center-of-mass reference frame. Solving 
(13.11) with the current value for the dark matter relic density Q^h"^ = 0.105 ± 0.008 [33.] 
gives: 

(^^12fXann) ^ (1-9 ± 0.2) X lO'^GcV-^, (3.2) 

for a range of dark matter masses between roughly lOGeV to lOOGeV and Xf between 



19.2 and 21.6, with about 0.4 thickness 34 



The value in (13. 2p for the dark matter annihilation cross-section translates into a 
relation between the parameters of a given theory entering the calculated expression of 
(fi2crann), heuce imposiug a constraint on these parameters which will limit the intervals 
of possible dark matter masses. This constraint can be exploited to examine aspects of 
the theory like perturbativity. For example, in our model, we can obtain via (13. 2p the 
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mutual coupling constant rj^i for given values of Aq , study its behavior as a function 
of rriQ and tell which dark-matter mass regions are consistent with perturbativity. Note 
that once the two mutual coupling constants Aq''^ and rj^f are perturbative, all the other 
physical coupling constants will be. In the study of this section, we choose the values 
Ag"^"* = 0.01 (very weak), 0.2 (weak) and 1 (large). We also let the two masses mo and 
nil stretch from O.lGeV to 120GeV, occasionally mo to 200GeV. Finally, note that we 
do not incorporate the uncertainty in (13.21) when imposing the relic-density constraint, 
something that is sufficient in view of the descriptive nature of this work. 



e = 10°, Ao'"^' =0.01, mi = lOGeV 




FIG. 1: 7]qi vs mo for small mi, small mixing and very small WIMP-Higgs coupling. 

The dark matter annihilation cross sections (times the relative speed) through all 
possible channels are given in the appendix. The quantity (f i2C"ann) is the sum of all these 
contributions. Imposing (fi20"ann) = 1.9 x 10~^GeV~^ dictates the behavior of 770^'*, which 
is displayed as a function of the dark matter mass mo. Of course, as the parameters are 
numerous, the behavior is bound to be rich and diverse. We cannot describe every bit of 
it. Also, one has to note from the outset that for a given set of values for the parameters, 
the solution to the relic-density constraint is not unique: besides positive real solutions 
(when they exist), we may find negative real or even complex solutions. It is beyond the 
scope of the present work to investigate the nature and behavior of all the solutions. We 
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= 10°, Ao*'*^ =0.01, mi = 30GeV 




FIG. 2: tJq^ vs mg for moderate mi, small mixing and very small WIMP-Higgs coupling. 

are only interested in finding the smallest positive real solution t]1^ when it exists, looking 
at its behavior and finding out when it is small enough to be perturbative. 

A. Small mixing angle and very weak dark matter — Higgs coupling 

Let us describe briefiy and only partly how the mutual Sq - Si coupling constant rj^f 
behaves as a function of the 5*0 mass mo- We start by a small mixing angle, say 6 = 10°, 
and a very weak mutual 5*0 - Higgs coupling constant, say Aq*^'' = 0.01. Let us also fix the 
5*1 mass first at the small value mi = lOGeV. The corresponding behavior of rjl^ versus 
mo is shown in Fig. [1] 

The range of mo shown is from O.lGeV to 200GeV, cut in four intervals to allow for 
'local' features to be displayed*. We see that the relic-density constraint on 5*0 annihilation 
has no positive real solution for mo < 1.3GeV, and so, with these very small masses, 5*0 
cannot be a dark matter candidate. In other words, for mi = lOGeV, the particle 5*0 
cannot annihilate into the lightest fermions only; inclusion of the c-quark is necessary. 
Note that right about mo — 1.3GeV, the c threshold, the mutual coupling constant 7701'' 



* A logplot in this descriptive study is not advisable. 
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0=lO°,V^^ = O.Ol,mo = O.2GeV 




FIG. 3: tJq^ vs mi for very light Si, small mixing and very small WIMP-Higgs coupling. 

starts at about 0.8, a value, while perturbative, that is roughly eighty-two-fold larger 
than the mutual 5*0 - Higgs coupling constant Aq'^''. Then r]^^ decreases, steeply first, 
more slowly as we cross the r mass towards the h mass. Just before mi/2, the coupling 
r]^i hops onto another solution branch that is just emerging from negative territory, gets 
back to the first one at precisely mi/2 as this latter carries now smaller values, and then 
jumps up again onto the second branch as the first crosses the mo-axis down. It goes 
up this branch with a moderate slope until mo becomes equal to mi, a value at which 
the 5*1 annihilation channel opens. Right beyond mi, there is a sudden fall to a value 
r^Q^^ ~ 0.0046 that is about half the value of \ and rj^i stays fiat till mo — 45GeV where 
it starts increasing, sharply after 60GeV. In the mass interval mo — 66GeV— 79GeV, there 
is a desert with no positive real solutions to the relic-density constraint, hence no viable 
dark matter candidate. Beyond mo — 79GeV, the mutual coupling constant iq^Q^ keeps 
increasing monotonously, with a small notch at the W mass and a less noticeable one 
at the Z mass. Note that for this value of mi (lOGeV), all values reached by r^Q^"* in 
the mass range considered, however large or small with respect to \ are perturbativily 
acceptable. 

Increasing mi to moderate values does not change the above qualitative features. As 
an illustration. Fig. |2] shows the behavior of function of mo for mi = 30GeV, 

keeping the mixing angle = 10°, still small, and the mutual 5*0 - Higgs coupling constant 
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G = 10°, Ao^^' = 0.01, mo = 1.4GeV(L) and 1.5GeV(R) 




FIG. 4: tJq^ versus mi for mo above r threshold. 

Ag'^'' = 0.01, still very weak. The first thing to note is that not all values of rj^^^ are 
perturbative. Indeed, tIq^ does not start until tuq ~ l.SGeV, but with the very large 
value"'' 89.8. It decreases very sharply right after, to 2.04 at about 1.6GeV. It continues 
to decrease with a pronounced change in the slope at the b threshold. Effects at the 
masses mi/2 and mi similar to those of figure [1] do occur here too. There is a desert 
that lies in this case in the mass interval 66.5GeV - 76.5GeV. At the upper bound, the 
coupling 7]^^ takes the value 2.15 and decreases very slowly till mo — 78.2GeV. Right 
after this mass, it plunges down to catch up with a solution branch that is just emerging 
from negative values. This solution branch increases steadily with two small notches at 
the W and Z masses. A similar global behavior occurs at other moderate mi masses, 
with varying local features. 

Because of the very-small-mo deserts described and visible on Fig. [H one may ask 
whether the model ever allows for very light dark matter. To look into this, we fix mo at 
small values and let mi vary. Take first mo = 0.2GeV and see Fig. |3l The allowed 5*0 
annihilation channels are the very light fermions e,u,d,fi and s, plus 5*1 when mi < rriQ. 



^ This feature is not displayed in figure [2] to avoid masking tlie other much smaller values taken by the 
mutual coupling. 
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^=10°,V^^ = 0.2,mi = 20GeV 




FIG. 5: tJq^ vs itiq for small mixing, moderate mi and WIMP-Higgs coupling. 

Note that we still have 9 = 10° and Ag'*^ = 0.01. Qualitatively, we notice that in fact, there 
are no solutions for mi < 0.2GeV (= mo here), a mass at which rj^^ takes the very small 
value ~ 0.003. It goes up a solution branch and leaves it at mi ~ 0.4GeV to descend 
on a second branch that enters negative territory at mi ~ 0.7GeV, forcing 7]^^ to return 
onto the first branch. There is an accelerated increase till mi ~ 5GeV, a value at which 
r^Qi^ ^0.5. And then a desert, no positive real solutions, no viable dark matter. 

Increasing mo until about 1.3GeV does not change these overall features: some 'move- 
ment' for very small values of mi and then an accelerated increase till reaching a desert 
with a lower bound that changes with m^. For example, the desert starts at mi ~ 6.8GeV 
for mo = 0.6GeV and mi ~ 7.3GeV for mo = 1.2GeV. Note that in all these cases where 
mo < 1.3GeV, all values of 7]^^ are perturbative. Therefore, the model can very well 
accommodate very light dark matter with a restricted range of 5*1 masses. 

However, the situation changes after the inclusion of the r annihilation channel. Indeed, 
as Fig. m shows, for mo = 1.4GeV, though the overall shape of the behavior of rj^^ as 
a function of mi is qualitatively the same, the desert threshold is pushed significantly 
higher, to mi ~ 20GeV. But more significantly, tiq^ starts to be larger than one already at 
mi ~ 17GeV, therefore loosing perturbativity. For mo = l.SGeV, the desert is effectively 
erased as we have a sudden jump to highly non-perturbative values of 7]^^ right after 
mi ~ 28GeV. Such a behavior stays with larger values of mo. But for mi < 20GeV (case 
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mo = l.SGeV), the values of t^qi ^ire smaller than one and physical use of the model is 
possible if needed. 



6'=10°,V^^ = 0.2,mi = 60GeV 




FIG. 6: tjq^ versus mo for heavy Si, small mixing and small WIMP-Higgs coupling. 



B. Small mixing angle and larger dark matter Higgs couplings 

What are the effects of the relic-density constraint when we vary the parameter A[f'*? 
Let us keep the Higgs - 5*1 mixing angle small {6 = 10°) and increase X^q\ first to 0.2 
and later to 1. For aJ,^^ = 0.2, Fi gure shows the behavior of function of the 

dark matter mass mo when mi = 20GeV. We see that rjl^ starts at mo — 1.4GeV with 
a value of about 1.95. It decreases with a sharp change of slope at the b threshold, then 
makes a sudden dive at about 5 GeV, a change of branch at mi/2 down till about 12GeV 
where it jumps up back onto the previous branch just before going to cross into negative 
territory. It drops sharply at mo = mi and then increases slowly until mo — 43.3GeV. 
Beyond, there is nothing, a desert. 

This is of course different from the situation of very small X'^^ like in Fig. [T]and Fig. |2] 
above: here we see some kind of natural dark-matter mass 'confinement' to small-moderate 
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viable-l- values. 

Still for Ag'^^ = 0.2 with mi = 60GeV this time, Fig. [H] shows rj^^ starting very high 
(~ 85GeV) at mg — l.SGeV, decreasing quickly with a first sudden drop at 2.7 GeV and 
a second one to zero at 26.2 GeV. A solution branch is then picked up - left briefiy at 
nil/ 2 ^ until 49 GeV and then nothing. What is peculiar here is that, in contrast with 
previous situations, the desert starts at a mass rriQ < mi, i.e., before the opening of the 
Si annihilation channel. In other words, the dark matter is annihilating into the light 
fermions only and the model is perturbatively viable in the range 20GeV - 49GeV. 



e - 10°, Ao^"^^ ^ 1, mi = 20GeV 




mo (GeV) 



FIG. 7: T^Qi versus mg for medium mi, small mixing and large WIMP-higgs coupling. 

The case A[,^^ = 1 with mi = 20GeV is displayed in Fig. [7l There are no solutions 
below mo — l.SGeV at which rjl^ ~ 1.80. From this value, tIq^^ slips down very quickly 
to pick up less abruptly when crossing the r threshold. There is a significant change in 
the slope at the crossing of the b mass. Note the absence of a solution at mi/2, which is 
a new feature, present for other values of mi not displayed here. Beyond mi/2, there is 
a slight change in the downward slope, a change of solution branch, and that goes until 
14.5GeV where t^qi'' jumps to catch up with the previous branch. It goes down this branch 
until about ISGeV where the desert starts. 

We have studied the behavior of 

^01 ^ 

function of mo for other values of mi between 
20GeV and lOOGeV while keeping 6 = 10° and = 1. For mi < 79.2GeV, the behavior 



Note that the values of 7]^^ for 1.6GeV < mg < 43.3GeV are all perturbative. 
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= 10°, A,/"*' = 1, mi = 79.3GeV 
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FIG. 8: tiq^ versus mo for heavy Si, small mixing and large WIMP-Higgs coupling. 



is qualitatively quite similar to that shown in Fig. [TJ but beyond this mass, there is a 
highly non-perturbative branch r]^^ jumps onto at small and moderate values of mo. This 
highly non-perturbative region stretches in size as mi increases. Fig. [H] displays this new 
feature. Note that on this figure, not all of the range of Vj^^ is shown in order to allow the 
small-coupling regions to be displayed; the high values of t]^^ are in the two thousands. 
Note also that it is the same highly non-perturbative solution branch r]^^ jumps onto for 
other large values of mi. 



C. Larger mixing angles 

Last in this descriptive study is to see the effects of larger values of the Si - Higgs 
mixing angle 9. We give it here the value 9 = 40° and tune back the mutual Sq - Higgs 
coupling constant X^^^ to the very small value 0.01. Figure |9] shows the behavior of tiq^ 
as a function of mo for mi = 20GeV. One recognizes features similar to those of the 
case 9 = 10°, though coming in different relative sizes. The very-small-mo desert ends at 
about 0.3GeV. There are by-now familiar features at the c and b masses, mi/2 and mi. 
Two relatively small forbidden intervals (deserts) appear for relatively large values of the 
dark matter mass: 67.3GeV — 70.9GeV and 79.4GeV — 90.8GeV. The W mass region is 
forbidden but there is action as we cross the Z mass. 

Other values of mi, not displayed because of space, behave similarly with an additional 
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= 40°, Ao*"** = 0.01, mi = 20GeV 




FIG. 9: iJq^ versus mo for moderate mi, moderate mixing and small WIMP-Higgs coupling. 

effect, namely, a sudden drop in slope at mg = {rrih + /«i)/2 coming from the ignition 
of 5*0 annihilation into Si and Higgs. We have also worked out the cases X^^^ = 0.2 
and 1 for 6* = 40°. The case X^^^ = 0.2 is displayed in Fig. dU] and presents differences 
with the corresponding small-mixing situation 6 = 10°. Indeed, for mi = 20GeV, the 
first feature we notice is a smoother behavior; compare with Fig. \5\ Here, 7]^^ starts at 
mo — 0.3GeV with the small value ~ 0.016 and goes up, faster at the c mass and with a 
small effect at the b mass. It increases very slowly until mi/2 and decreases very slowly 
until mo = mi, and then there is a sudden change of branch followed immediately by a 
desert^. So here too the model naturally confines the mass of a viable dark matter to 
small-moderate values, a dark matter particle annihilating into light fermions only. What 
is also noticeable is that there is stability of rj^f around the value of X'^^ in the interval 
l.SGeV - 20GeV (= mi here). 

The case mi = 60GeV presents also overall similarities as well as noticeable differences 
with the corresponding case 6 = 10°, see Fig. |6l The first difference is that all values of ?7oi'* 
are perturbative. This latter starts at mo — 1.4GeV with the value ~ 0.75, goes down and 
jumps to catch up with another solution branch emerging from negative territory when 



§ Except for the very tiny interval 78.5GeV - 79.0GeV not displayed on Fig. [TOj 
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Q = 40°, Ao^^^ = 0.2, mi = 20GeV(L) and 60GeV(R) 




FIG. 10: rjQ^ versus uiq for moderate (L) and large (R) mi, large mixing and moderate WIMP- 
Higgs coupling. 

crossing the r mass. It increases, kicking up when crossing the 6-quark mass. It changes 
slope down at mi/2 and goes to zero at about 51GeV. It jumps up onto another branch 
that goes down to zero also at about 58.6GeV, just below mi, and then there is a desert, 
except for the small interval 76.3GeV — 80.5GeV. 

The case Aq^^ = 1 is shown in Fig. [TTJ Global similarities with the previous case are 
apparent. All values of Vj^^ are perturbative and the mass range is naturally confined to 
the interval 0.2GeV - 20GeV for mi = 20GeV, and 1.4GeV-52.3GeV for mi = 60GeV. 
We note action at the usual masses and, in particular, we see there are no solutions at 
mo = mi/2 like in the case 9 = 10°. We note here too the quasi- constancy of ?7ot^ for most 
of the available range. 

Finally, we note that we have worked out larger mixing angles, notably 6 = 75°. In 
general, these cases do not display any new features worth discussing: the overall behavior 
is mostly similar to what we have seen, with expected relative variations in size. 
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9 = 40°, Ao*^' = 1, mi = 20GeV(L) and 60GeV(R) 
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FIG. 11: r]l^ versus ttlq for moderate (L) and large (R) mi, large mixing and large WIMP-Higgs 
coupling. 



IV. DARK-MATTER DIRECT DETECTION 

Experiments like CDMS II |27|? ], XENON 10/100 (sl, DAMA/LIBRA Q and 
CoGeNT jsj search directly for a dark matter signal. Such a signal would typically come 
from the elastic scattering of a dark matter WIMP off a non-relativistic nucleon target. 
However, throughout the years, such experiments have not yet detected an unambiguous 
signal, but rather yielded increasingly stringent exclusion bounds on the dark matter - 
nucleon elastic-scattering total cross section a^et in terms of the dark matter mass mo- 

In order for a theoretical dark-matter model to be viable, it has to satisfy these bounds. 
It is therefore natural to inquire whether the model we present in this work has any 
capacity of describing dark matter. Hence, we have to calculate cxdet as a function of 
mo for different values of the parameters {6,X'^\mi) and project its behavior against 
the experimental bounds. We will limit ourselves to the region O.lGeV - lOOGeV as 
we are interested in light dark matter. As experimental bounds, we will use the results 
from CDMSII and XENONIOO, as well as the future projections of SuperCDMS Q and 
XENONItH. The results of CoGeNT, DAMA/LIBRA and CRESST will be discussed 
elsewhere j35|. As the figures below show jsGj, in the region of our interest, XENONIOO 
is only slightly tighter than CDMSII, SuperCDMS significantly lower and XENONIT the 
most stringent by far. But it is important to note that all these results loose reasonable 
predictability in the very light sector, say below 5GeV. 

The scattering of 5*0 off a SM fermion / occurs via the t-channel exchange of the 
SM higgs and 5"!. In the non-relativistic limit, the effective Lagrangian describing this 
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interaction reads 



afffSo, 



(4.1) 



where 



2v 



Aq^^ cos 6 7]^^' sin 6 



(3), 



mt 



In this case the total cross section for this process is given by 



m 



f 



An {rrif + mo) v"^ 



Xf^ cos 6 tJq^^' sin 6 



(3) 



n 2 



mt 



(4.2) 



(4.3) 



At the nucleon level, the effective interaction between a nucleon N = p or n and 5*0 has 
the form 



4o-N = aNNNSl 
where the effective nucleon-^o coupling constants is given by: 



{mN - 



Aq^'* cos 6 riQi' sin 6 



(3), 



mt 



2 



(4.4) 



(4.5) 



In this relation, m^r is the nucleon mass and the baryon mass in the chiral limit 25 . 
The total cross section for non-relativistic Sq - N elastic scattering is therefore: 



C"det = 0■SoN-^SoN 



m% [mN - 1^3)' 
An [m^ + itiq)'^ f ^ 



X^o^ cose 
ml 



(3) • a 
tIqi sm u 



mt 



(4.6) 



The rest of this section is devoted to a brief discussion of the behavior of cxdet as 
a function of mo. We will of course impose the relic-density constraint on the dark 
matter annihilation cross section (13. 2p . But in addition, we will require that the coupling 
constants are perturbative, and so impose the additional requirement < rj^^ < 1. Also, 
here too, the choices of the sets of values of the parameters {6, Xl^\mi) can by no means be 
exhaustive but only indicative. Furthermore, though a detailed description of the behavior 
of (Tdet could be interesting in its own right, we will refrain from doing so in this work as 
there is no need for it, and content ourselves with mentioning overall features and trends. 
Generally, as m^ increases, the detection cross section adot starts from high values, slopes 
down to minima that depend on the parameters and then picks up moderately. There are 
features and action at the usual mass thresholds, with varying sizes and shapes. Excluded 
regions are there, those coming from the relic-density constraint and new ones originating 
from the additional perturbativity requirement. Close to the upper boundary of the mass 
interval considered in this study, there is no universal behavior to mention as in some cases 



19 




FIG. 12: Elastic N — Sq scattering cross-section as a function of mo for moderate mi, small 
mixing and small WIMP-Higgs coupling. 

(Jdet will increase monotonously and, in some others, it will decrease or 'not be there' at 
all. Let us finally remark that the logplots below may not show these general features 
clearly as these latter are generally distorted. 

Let us start with the small Higgs - 5*1 mixing angle 6 = 10° and the very weak 
mutual 5*0 - Higgs coupling Xl^^ = 0.01. Fig. [12] shows the behavior of o"det versus mo 
in the case mi = 20GeV. We see that for the two mass intervals 20GeV — 65GeV and 
75GeV — lOOGeV, plus an almost singled-out peaks at uiq = mi/2, the elastic scattering 
cross section is below the projected sensitivity of SuperCDMS. However, XENONIT will 
probe all the these masses , except mo — 58 GeV and 85 GeV. 

Also, as we see in Fig. [121 most of the mass range for very light dark matter is excluded 
for these values of the parameters. Is this systematic? In general, smaller values of mi 
drive the predictability ranges to the lighter sector of the dark matter masses. Figure 
[T3] illustrates this pattern. We have taken mi = 5GeV, just above the lighter quarks 
threshold. In the small-mass region, we see that SuperCDMS is passed in the range 
5GeV — 30GeV. Again , all this mass ranges will be probed by XENONIT experiment, 
except a sharp peak at mo = mi/2 = 2.5GeV, but for such a very light mass, the 
experimental results are not without ambiguity. 

Reversely, increasing mi shuts down possibilities for very light dark matter and thins 
the intervals as it drives the predicted masses to larger values. For instance, in figure [Ml 
where mi = 40GeV, in addition to the peak at mi/2 that crosses SuperCDMS but not 
XENONIT, we see acceptable masses in the ranges 40GeV - 65GeV and 78GeV - up. 
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e = 5\ Ao^"^^ = 0.01, mi = 5GeV 
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FIG. 13: Elastic N — Sq scattering cross-section as a function of niQ for light Si, small mixing 
and small WIMP-Higgs coupling. 

Here too the intervals narrow as we descend, surviving XENONIT as spiked peaks at 62 
GeV and around 95 GeV. 

A larger mutual coupling constant Ag^'' has the general effect of squeezing the acceptable 
intervals of uiq by pushing the values of Udet up. As an illustration, see figure [15] where 
we have taken Ag^'' = 0.2 and a larger value of mi = 60GeV. In this example, already 
XENONIOO excludes all the masses below 100 GeV except a relatively narrow peak at 
mi/2. 

Increasing the mixing angle 6 has also the general effect of increasing the value of 
(Tdet- Figure [16] shows this trend for 6 = 40°; compare with Fig. [121 The only allowed 
masses by the current bounds of CDMSII and XENONIOO are between 20 GeV and 50 
GeV, the narrow interval around mi/2, and another very sharp one, at about 94GeV. The 
projected sensitivity of XENONIT will probe all mass ranges except those at mg ~ 30GeV 
and 94GeV. 

Finally, it happens that there are regions of the parameters for which the model has no 
predictability. See figure [17] for illustration. We have combined the effects of increasing 
the values of the two parameters Ag'^'' and mi. As we see, we barely get something at 
mi/2 that cannot even cross XENONIOO down to SuperCDMS. 
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e = 10°, Ao^"^^ = 0.01, mi = 40GeV 
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FIG. 14: Elastic N — Sq scattering cross-section as a function of niQ for medium mi, small 
mixing and small WIMP-Higgs coupling. 



V. CONCLUDING REMARKS 

In this work, we presented a plausible scenario for light cold dark matter, (for masses 
lighter than 100 GeV). This latter consists in enlarging the Standard Model with two 
gauge-singlet Z2-symmetric scalar fields. One is the dark matter field 5*0, stable, while 
the other undergoes spontaneous symmetry breaking, resulting in the physical field 5*1. 
This opens additional channels through which 5*0 can annihilate, hence a reducing its 
number density. The model is parametrized by three quantities: the physical mutual 
coupling constant Aq"^^ between 5*0 and the Higgs, the mixing angle 6 between 5*1 and the 
Higgs and the mass mi of the particle 5*1. We first imposed on 5*0 annihilation cross 
section the constraint from the observed dark-matter relic density and studied its effects 
through the behavior of the physical mutual coupling constant r/oi^ between 5*0 and 5*1 as 
a function of the dark matter mass rriQ. Apart from forbidden regions (deserts) or others 
where perturbativity is lost, we find that for most values of the three parameters, there 
are viable solutions in the small-moderate masses of the dark matter sector. Deserts are 
found for most of the ranges of the parameters whereas perturbativity is lost mainly for 
larger values of mi. Through the behavior of rj^^ , we could see the mass thresholds which 
mostly affect the annihilation of dark matter, and these are at the c, r and b masses, as 
well as mi/2 and mi. 

The current experimental bounds from CDMSII and XENONIOO put a strong con- 
straint on the 5*0 masses in the range between 10 to 20 GeV. For small values of mi, very 
light dark matter is viable, with a mass as small as one GeV. This is of course useful for 
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FIG. 15: Elastic N — Sq scattering cross-section as a function of tuq for heavy Si, small mixing 
and moderate WIMP-Higgs coupling. 
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understanding the results of the experiments DAMA/LIBRA, CoGeNT , CRESST 
as well as the recent data of the Fermi Gamma Ray Space Telescope [3] • The projected 
sensitivity of future WIMP direct searches such as XENONIT will probe all the 5*0 masses 
between 5 GeV and 100 GeV. 

The next step to take is to test the model against the phenomenological constraints. 
Indeed, one important feature of the model is that it mixes the 5*1 field with the Higgs. 
This must have implications on the Higgs detection through the measurable channels. 
Current experimental bounds from LEP II data can be used to constrain our mixing angle 
9, and possibly other parameters. In addition, a very light 5*0 and/or will contribute to 
the invisible decay of J/ ip and T mesons and can lead to a significant branching fraction. 
These constraints can be injected back into the model and restrain further its domain of 



validity. These issues are under current investigation [35 . 

Also, in this work, the 5*1 vacuum expectation value vi was taken equal to lOOGeV, 
but a priori, nothing prevents us from considering other scales. However, taking vi much 
larger than the electro-weak scale requires 7]"^^ to be very tiny , which will result in the 
suppression of the crucial annihilation channel SqSq — )■ SiSi . Also, we have fixed the 
Higgs mass to = 138GeV, which is consistent with the current acceptable experimental 



bounds [30|. Yet, it can be useful to ask here too what the effect of changing this mass 
would be. 

Finally, in this study, besides the dark matter field 5*0, only one extra field has been 
considered. Naturally, one can generalize the investigation to include such fields and 
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FIG. 16: Elastic N — Sq scattering cross-section as a function of for moderate mi, large 
mixing and small WIMP-Higgs coupling. 

discuss the cosmology and particle phenomenology in terms of A^. It just happens that 
the model is rich enough to open new possibilities in the quest of dark matter worth 
pursuing. 



Appendix A: Dark matter annihilation cross-sections 

The cross-sections related to the annihilation 5*0 into the scalar particles are as follows. 
For the hh channel, we have: 
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FIG. 17: elastic cross-section Uei as a function of Sq mass for heavy Si, small mixing and 
relatively large WIMP-Higgs coupling. 



The function is the step function. For the 5*1 5*1 channel, we have the result: 
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For the hSi channel, we have: 
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The annihilation cross-section into fermions is: 
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Last, the annihilation cross-section into Z's is: 
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The quantities = Jti^F/j and ei = miFi are regulators at the respective resonances. 
The decay rates and Fi are calculable in perturbation theory. We have for h: 
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For 5*1, we have similar expressions: 
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where Nc is equal to 1 for leptons and 3 for quarks. 
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